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1 Introduction 

The estimation of the probability of an 
individual nearby lightning stroke 
contacting the surface within a specified 
distance of a specified spaceport 
processing facility is important to allow 
engineers to decide if inspection of 
electronics for damage from induced 
currents from that stroke is warranted. If 
induced current damage has occurred, 
inspections of the electronics are critical 
to identify required fixes and avoid 
degraded performance or failure of the 
satellite or space launch vehicle. 
However, inspections are costly both 
financially and in terms of delayed 
processing for space launch activities. As 
such, it is important these inspections be 
avoided if not needed. At Kennedy Space 
Center (KSC) and Cape Canaveral Air 
Force Station (CCAFS), one of the main 
purposes of the Four Dimensional 
Lightning Surveillance System (4DLSS) 
[8], [9] is detection of nearby strokes and 
determination of their peak current to 
support those inspection decisions 
[5], [6], [10]. The high frequency of 
occurrence of lightning in East Central 
Florida combined with the large amount of 
complex sensitive electronics in satellite 
payloads, space launch vehicles, and 
associated facilities make those decisions 
critically important to space launch 
processing. While 4DLSS provides the 


data for 95 th or 99th percentile location 
error ellipses for the best location for each 
stroke, depending on customer, it has not 
been able to provide the probability for the 
stroke being within a customer specified 
distance of a point of interest. This paper 
presents a new method to convert the 
4DLSS error of 50th percentile location 
error ellipse for best location of any stroke 
into the probability that the stroke was 
within any radius of any facility at 
CCAFS/KSC. This new facility-centric 
technique is a significant improvement 
over the stroke-centric location error 
ellipses the 45 th Weather Squadron (45 
WS) was providing. This technique is 
adapted from a method of calculating the 
probability of debris collision with 
spacecraft. [1], [2], [3], [6] A table of 

abbreviations used in this paper is 
available in Section 7, Table 2. 

2 Background 

In spacecraft collision probability and 
the other applications, at the instant of 
“nominal” closest approach, the position 
uncertainty of the collision object relative 
to the asset is described by a bivariate 
Gaussian probability density function 
(pdf).[1],[2],[3] 
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where a x and o z = the standard deviations 
of x and z, pxz = correlation coefficient of x 
and z, x and z are the designations for the 
rectangular coordinates in the collision 
plane. 

The probability of collision is given by 
the two-dimensional integral, where A is 
the collision cross-sectional area which is 
a circle with radius, r A .[3] 


P = ^f 2 (x,z)dxdz 

A 

There is no known analytical solution to 
the above integral when the two standard 


deviations o x and o z are not equal. The 
solution is based either on transforming 
the two dimensional Gaussian probability 
distribution function (pdf) to a one- 
dimensional Rician pdf and using the 
concept of equivalent areas or by 
performing a numerical integration of the 
two dimensional Gaussian pdf. [1], [2], [3] 

The geometry used for spaceflight 
collision probability can also be used for 
estimation of the probability of an 
individual nearby lightning stroke 
contacting the surface within a specified 
distance of a specified point of interest as 
shown in Figure 1 below. Both solution 
methods, numerical integration as well as 
the analytical method of equivalent areas 
using the Rician pdf, will be analyzed in 
the next section. 



Figure 1. Schematic diagram of the angles used in probability calculation for a sample 
lightning location error ellipse, a is the heading of the semi-major axis of the lightning 
location uncertainty ellipse from true north. 0 is the angle between the semi-major axis 
of the lightning location uncertainty ellipse and line connecting the center of the lightning 
uncertainty ellipse and the center of the area of interest. 


3 Evaluation 


3. 1 Test Set 1 


The probability that any lightning strike 
is within any radius of any point of interest 
would be extremely difficult to estimate 
intuitively. As a result, given the high 
impact resultant decisions, the tool must 
be extensively tested. Three major types 
of tests were conducted and are 
discussed in the following sections: 
1) known mathematical solutions (Test 
Set 1 ), 2) expected behavior as single 
parameters are varied (Test Set 2), and 
3) examination of real-world events (Test 
Set 3 and Test Set 4). The new technique 
passed all of the tests. 


The first set of testing compared the 
probability calculated by the program to 
the corresponding circular probability from 
the CRC Handbook of Tables for 
Probability and Statistics. [4] Table 1 
shows, for various inputs, the calculated 
probability and the CRC Handbook 
probability. The values matched to within 
a tenth of a percent. These differences in 
the final digit may be due to round-off 
error. 
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Table 1. Calculated probability vs. CRC Handbook probability for various inputs. The 
abbreviations are defined as follows. 


Abbreviation 

Definition 

conf 

confidence 

CE 

confidence ellipse 

lat 

latitude 

Ion 

longitude 

nmi 

nautical miles 

POI 

point of interest 






















































3.2 Test Set 2 

The second type of testing involved 
plotting the calculated probabilities as 
particular inputs were varied while holding 
the other inputs constant and comparing 
them to an independently coded program 
written by Dr. F. Kenneth Chan of the 
Aerospace Corporation and the author of 
"Spacecraft Collision Probability". [3] Also 
tested was the difference between using a 
numerical integration technique for 
calculating the probability versus an 
analytical technique, shown as “Rician” in 
the results below. The analytical or Rician 
technique involves transforming the two 
dimensional Gaussian pdf to a one- 
dimensional Rician pdf and using the 
concept of equivalent areas to calculate 
the probability.[1],[2],[3] The results are 
shown in Figures 2 through 5 and 8 
below. The data used to generate these 
figures are in Table 2. Note that results 
using the 45 th Weather Squadron (45WS) 
and Chan’s program match almost exactly 
regardless of integration method used. 
Probability calculations are much faster 
using the analytical technique as opposed 
to the numerical integration technique. 
Therefore, it was of interest to understand 
the conditions in which the analytical 
technique performed well compared to the 
numerical integration technique. The 
numerical integration technique and the 
analytical integration technique tend to 
diverge as the ratio of the semi-major axis 
to semi-minor axis increases and also as 
the orientation angle of the ellipse 
approaches 0 or 180 degrees. 

Figure 2 shows the change in 
probability as a result of changing the 
radius around the point of interest while 
holding all other parameters constant. 
Chan's probability calculated using both 
the numerical integration technique as 
well as the analytical (Rician) technique is 
compared to the probability calculated 


using the 45WS program. The worst case 
probability difference between methods is 
0.25% at a radius of 2 nautical miles 
around the point of interest. Chan’s 
probability using both techniques matches 
the 45WS probability exactly at all radii. 

Figure 3 shows the change in 
probability as a result of changing the 
latitude of the strike from the point of 
interest while holding all other parameters 
constant. Chan's probability calculated 
using both the numerical integration 
technique as well as the analytical 
(Rician) technique is compared to the 
probability calculated using the 45WS 
program. The probability follows a 
Gaussian curve and reaches a maximum 
when the uncertainty ellipse is at its 
closest point of approach to the point of 
interest, as expected. The worst case 
probability difference between methods is 
0.06 where the lightning stroke is at 
latitude of 28.61 62°N, which is about 0.5 
nautical miles away from the point of 
interest. Chan’s probability using both 
techniques matches the 45WS probability 
exactly at all latitudes. 

Figure 4 shows the change in 
probability as a result of changing the 
longitude of the strike from the point of 
interest while holding all other parameters 
constant. Chan's probability calculated 
using both the numerical integration 
technique as well as the analytical 
(Rician) technique is compared to the 
probability calculated using the 45WS 
program. The probabilities follow a 
Gaussian curve and reach a maximum 
when the uncertainty ellipse is at its 
closest point of approach to the point of 
interest, as expected. The worst case 
probability difference between methods is 
0.08 where the lightning stroke is at a 
longitude of 80.5961 'W, which is about 
0.4 nautical miles away from the point of 
interest. Chan’s probability using both 


techniques matches the 45WS probability 
exactly at all longitudes. 

Figure 5 shows the change in 
probability as a result of changing the 
heading from true north of the semi-major 
axis of the lightning uncertainty ellipse 
while holding all other parameters 
constant. Chan’s probability calculated 
using both the numerical integration 
technique as well as the analytical 
(Rician) technique is compared to the 
probability calculated using the 45WS 
program. The center of the stroke 
uncertainty ellipse is located about 0.5 
nautical miles away from the point of 
interest. The probabilities show a roughly 
sinusoidal pattern as more then less then 
more again of the ellipse rotates into, out 
of, then into the area around the point of 
interest. However the difference in 
probability between the two integration 
techniques is enhanced as the ellipse is 
rotated. The worst case probability 
difference between methods is 0.28 
where the lightning stroke heading is at an 
angle of 0° or 180° from true north. 

Figures 6 and 7 show a Google Maps 
visualization of the 99% confidence 
lightning uncertainty ellipse as it is rotated 
from 180° to 90° from true north. The 
lightning uncertainty ellipse at a heading 
of 90° from true north is the rotation angle 
at which there is no difference in the 
probability calculated by the numerical 
integration technique and the analytical 
(Rician) technique. Chan’s probability 


using both techniques matches the 45WS 
probability exactly at all angles. 

Figure 8 shows the change in 
probability as a result of varying the 
aspect ratio (length of semi-major 
axis/length of semi-minor axis) of the 
lightning uncertainty ellipse from 1.5 to 11 
with the strike point close to the point of 
interest while holding all other parameters 
constant. Chan's probability calculated 
using both the numerical integration 
technique as well as the analytical 
(Rician) technique is compared to the 
probability calculated using the 45WS 
program. The probability becomes less as 
the aspect ratio of the uncertainty ellipse 
is larger. However the difference in 
probability between the two integration 
techniques is enhanced as the aspect 
ratio is increased. The worst case 
probability difference between methods is 
0.07 where the aspect ratio is 8. Chan’s 
probability using both techniques matches 
the 45WS probability exactly at all aspect 
ratios. 

In light of the results of the differences 
between calculations (numerical 
integration vs. analytical (Rician) method, 
the 45WS has decided to use the 
numerical integration technique to 
calculate probabilities. Although the 
program run time is longer using the 
numerical integration technique, the 
accuracy improvements justify the longer/ 
calculation time. 
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Table 2. Input values used for scenarios shown in Figures 2 through 8. The 
abbreviations are the same as defined in Table 1. 
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igure 2. Change in probability as a 
result of changing the POI radius while 
holding all other parameters constant. 
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Figure 3. Change in probability as a 
result of changing the strike lat from the 
POI while holding all other parameters 
constant. 
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r igure 4. Change in probability as a 
result of changing the strike Ion of the 
strike from the POI while holding all other 
parameters constant. 




result of changing the semi-major axis 
heading of the lightning uncertainty ellipse 
while holding all other parameters 
constant. 



Figure 6. Google Maps visualization of a 
lightning uncertainty ellipse overlaid on 
the radius around the POI with a semi- 
major axis heading of 180° as graphed in 
Figure 5. 



Figure 7. Google Maps visualization of a 


lightning uncertainty ellipse overlaid on 
the radius around the POI with a semi- 
major axis heading of 90° as graphed in 
Figure 5. 
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Figure 8. Change in probability as a 
result of varying the aspect ratio (length of 
semi-major axis/length of semi-minor axis) 
of the lightning uncertainty ellipse from 1.5 
to 1 1 with the strike point close to the POI 
while holding all other parameters 
constant. 


3.3 Test Set 3 

The third type of testing analyzed six 
real-world lightning strikes near Space 
Launch Complex 39A on 3 August 2009. 
Figure 9 shows the spreadsheet used to 
generate the lightning report for those six 
strikes. Additional data on these strikes 
are in Table 3. These strikes were 
selected because the closest point on the 
lightning position uncertainty ellipse was 
within 0.45 nautical miles of Launch 
Complex 39A. Figures 10 through 15 are 
Google Maps depictions of these six 
strokes. The probabilities for a small area 
around a facility, even for a nearby stroke, 
may appear to be surprisingly low. For 
example, one strike just 0.65 nautical 
miles away (Figure 14) had only a 0.7% 
probability of being within the 
0.45 nautical mile radius of Launch 
Complex 39A. All calculated probabilities 
are consistent with these real-world 
events. 
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Figure 9. Sample of lightning strikes where the closest point on the lightning position 
uncertainty ellipse was within 0.45 nmi of Launch Complex 39A on 3 August 2009. 
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Table 2. Input values used for scenarios shown in Figures 10 through 15. The 
abbreviations are the same as defined in Table 1 . 



Figure 10. Google Maps visualization of 
the 99% conf uncertainty ellipse for the 
closest lightning strikes to Complex 39A 
on 03 August 2009. There is a 45.9% 
probability that the lightning strike 
occurred within the 0.45 nmi radius. 



Figure 11. Google Maps visualization of 
the 99% conf uncertainty ellipse for one of 
the closest lightning strikes to Complex 
39A on 03 August 2009. The center of 
the ellipse was within the 0.45 nmi radius. 
There is a 54.4% probability that the 
lightning occurred within that radius. 





Figure 12. Google Maps visualization of 
the 99% conf uncertainty ellipse for a 
lightning strike near Complex 39A on 03 
August 2009. Figure 12 shows a 
probability of 10.4% of the lightning strike 
occurring within the 0.45 nmi radius. 



Figure 13. Google Maps visualization of 
the 99% conf uncertainty ellipse for a 
lightning strike near Complex 39A on 03 
August 2009. Figure 13 shows a 
probability of 6.4% of the lightning strike 
occurring within the 0.45 nmi radius. 



Figure 14. Google Maps visualization of 
the 99% confidence uncertainty ellipse for 
nearby lightning strike to Complex 39A on 
03 August 2009. Figure 14 shows a 
probability of 0.7% of the lightning strike 
occurring within the 0.45 nmi radius. 



Figure 15. Google Maps visualization of 
the 99% conf uncertainty ellipse for the a 
lightning strike near Complex 39A on 03 
August 2009. Figure 15 shows a 
probability of 7.3% of the lightning strike 
occurring within the 0.45 nmi radius. 



3.4 Test Set 4 

The fourth type of testing analyzed six 
additional real-world lightning strikes in 
and around Space Launch Complexes 
39A and 39B. These examples were 
generated based on an Electromagnetic 
Environmental Effects (EEE) panel 
request to run several case studies where 
there was camera verification of a 
lightning strike in the vicinity of Launch 
Complex 39A or 39B. The data used for 
this analysis are in Table 4. Both CGLSS 
and NLDN cases were examined, 
depending upon which sensor system 
recorded the stroke. Figures 16 through 
21 show the probability results from these 
cases. As in section 3.3, all calculated 
probabilities were consistent with these 
additional real-world events. 



Figure 16. Illustrates a probability of 
92.1% of a lightning strike of amplitude - 
38.9 kA detected by CGLSS occurring 
0.32 nmi from the center of Launch 
Complex 39A on 8/16/2009. 
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Table 4. Input values used for scenarios shown in Figures 16 through 2 


. The 


abbreviations are the same as defined in Table 1 . 



Figure 17. Illustrates a probability of 
72.1% of a lightning strike of amplitude 
-43.0 kA detected by NLDN occurring 
0.26 nmi from the center of Launch 
Complex 39A on 8/16/2009. 



Figure 19. Illustrates a probability of 
97.2% of a lightning strike of amplitude - 
39.5 kA detected by CGLSS occurring 
0.12 nmi from the center of Launch 
Complex 39A on 7/21/2008. 



Figure 18. Illustrates a probability of 
77.7% of a lightning strike of amplitude 
-71.4 kA detected by NLDN occurring 
0.28 nautical miles from the center of 
Launch Complex 39A on 10/14/2009. 



Figure 20. Illustrates a probability of 
99.999975% of a lightning strike of 
amplitude -18.9 kA detected by CGLSS 
occurring 0.03 nmi from the center of 
Launch Complex 39B on 6/27/2009. 





Figure 21. Illustrates a probability of 
99.999925% of a lightning strike of 
amplitude -21.7 kA detected by CGLSS 
occurring 0.04 nmi from the center of 
Launch Complex 39B on 6/27/2009. 

4 Summary 

A technique has been developed to 
calculate the probability that any nearby 
lightning stroke is within any radius of any 
point of interest. In practice, this provides 
the probability that a nearby lightning 
stroke was within a key distance of a 
facility, rather than the error ellipses 
centered on the stroke. This process 
takes the current bivariate Gaussian 
distribution of probability density provided 
by the current lightning location error 
ellipse for the most likely location of a 
lightning stroke and integrates it to get the 
probability that the stroke is inside any 
specified radius. This new facility-centric 
technique will be much more useful to the 
space launch customers and may 
supersede the lightning error ellipse 
approach discussed in [5], [6]. 
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